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TECHNIQUES FOR EXAMINING 


STATISTICAL AND POWER- SPECTRAL PROPERTIES 

OF RANDOM TIME HISTORIES* 

By Herbert A. Leybold 
Langley Research Center 


SUMMARY 


A technique is presented for digitally generating random time histories 
having any desired shaped power spectra. Four random time histories having 
different statistical and power- spectral properties have been generated and 
analyzed to determine their instantaneous mean and amplitude distributions. In 
each, the distribution of instantaneous means could be approximated by a normal 
or Gaussian distribution and the distribution of instantaneous amplitudes could 
be approximated by the sum of a Rayleigh distribution and a normal distribution. 
An attempt was made to relate the coefficients of the equations used to repre- 
sent the distributions of means and amplitudes to the power- spectral properties 
of the generated time histories. Two of the coefficients could be related to 
the power- spectral properties of the time histories. The remaining two coeffi- 
cients were empirically determined since no apparent relationship was found 
between these coefficients and the power- spectral properties of the generated 
random time histories. 


INTRODUCTION 


Many of the loads encountered by aircraft and missiles are random in 
nature and, consequently, are usually described statistically. In order to 
reduce the mathematical complexity in utilizing such a description in analyzing 
the response of structures to loads, most investigators have made simplifying 
assumptions about the statistics of the random-load history (ref. l). As an 
example, in fatigue studies the statistics of the load peaks are usually used. 
These statistics are obtainable either by actually counting the peak loads at 
various levels or by a relationship developed by Rice (ref. 2 ) which relates 
the peak load distribution to the power spectrum of the random load-time history. 
When programing fatigue tests, all peak loads are usually applied about a com- 
mon mean load. In general, this mean load is representative of the overall mean 

*The information presented herein was offered as a thesis in partial ful- 
fillment of the requirements for the degree of Master of Science in Engineering 
Mechanics, Virginia Polytechnic Institute, Blacksburg, Virginia, May 1965* 


of the random load-time history from which the peak load distribution was 
derived. A variation in mean load can have an effect on fatigue life (ref. 3). 
Thus, it appears that a statistical description of both the distribution of 
instantaneous mean loads (i.e., average of two successive peak loads) and asso- 
ciated instantaneous amplitude distributions (i.e., difference between peak 
load and instantaneous mean load) would be more useful than the peak load dis- 
tribution alone for studying fatigue under random loading. The instantaneous 
mpan load distributions and associated instantaneous amplitude distributions are 
hereafter referred to as simply the mean and amplitude distributions. These 
distributions can be obtained by actually counting the instantaneous means and 
instantaneous amplitudes, but there is no known relationship between these dis- 
tributions and the power spectrum as was the case for the peak load 
distribution. 

In the present investigation, an attempt is made to develop an empirical 
relationship between the power- spectral properties of a given random time his- 
tory and the mean and amplitude distributions of this time history. This was 
done by digitally generating four random time histories with different power- 
spectral properties and counting the means and amplitudes in order to determine 
their distributions for each of the time histories generated. Equations 
describing the mean and amplitude distributions are developed and an attempt is 
made to relate the coefficients of these equations to the power- spectral prop- 
erties of the generated time histories. 


SYMBOLS 


A amplitude of periodic function of time 

aft filter factors or Fourier coefficients 

frequency-response function 
f frequency, cps 

fp Nyquist or folding frequency, cps 

f a number of amplitudes counted in interval yp -g y ^ y 1+1 

f^ number of amplitudes counted which exceed y = y^ 

f c computed frequency of occurrence of amplitudes in interval 

^i = y = yi+i 

f Q number of times per second zero axis is crossed with positive slope 

f number of positive peaks per second 

Jr 
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fy number of times per second that value y is exceeded 

fy^ computed number of amplitudes which exceed y = 

Kp,K g constants of peak probability distributions 

raw spectral density estimate 

1% number of amplitudes normally distributed about specified mean value 

% number of amplitudes distributed according to Rayleigh distribution 

about specified mean value 

Np number of positive amplitudes about specified mean value 

Pj^ normal probability 

Pir Rayleigh probability 

Pp(y) probability that peak will exceed given value of y 
p modified random number 

R(t) covariance function or autocorrelation function of continuous variable 

R]y generated random number 

Rp covariance or autocorrelation of discrete set of values 

t time, sec 

At uniform interval of time, sec 

X normal deviate 

Y(t) filtered time history 

Yi discrete set of values obtained by sampling filtered time history 

Y(t) at uniform intervals of time At 

y(t) original time history 

y^ discrete set of values obtained by sampling original time history 

y(t) at uniform intervals of time At 

z standard variable, y/Oy 

a, p dummy variables 
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°N 


°R 


°y 


2 

a. 

V 

a.? 

y 


coefficient of normal probability- 
coefficient of Rayleigh probability 

standard deviation or root-mean- square (rms) value of y(t) 
mean square value of y(t) 

mean square value of first derivative of y(t) 
mean square value of second derivative of y(t) 


0(f) power spectral density of a continuous variable 

0k smoothed spectral density estimate 

Matrix notations: 

[] square matrix 



column matrix 

A dot over a variable indicates differentiation with respect to time. 
A bar over a term indicates the mean value of the term. 


GENERATION OF RANDOM TIME HISTORIES 


A digital random time history having the properties of band-limited white 
noise was generated and used as the input to several linear systems, each 
having significantly different frequency- response characteristics. The output 
responses obtained were used to calculate power- spectral properties and also to 
obtain the distributions of means and amplitudes. 

A brief outline of the procedures used to simulate digitally random time 
histories having different shaped power spectra is as follows. A more detailed 
discussion of each of the following steps will appear in subsequent paragraphs. 

1. Random numbers having a uniform probability distribution were generated. 


2. The uniformly distributed random numbers were then transformed into a 
normal or Gaussian distribution having a mean of zero and a variance of one* 

The numbers obtained were assumed to be samples taken at l/2-sec intervals from 
a continuous record- A power spectrum was calculated by using these numbers 
and was found to be essentially flat. The normally distributed numbers will be 
used as the input to a linear system. 

3- In order to determine the frequency response of the linear system the 
following equation was used: 


0in 



^out 


where is the power spectrum of the normally distributed numbers obtained 
in step 2 above and 0 ou t is the desired shaped power spectrum. Knowing both 
the input and the desired output power spectrums, the magnitude of the fre- 


quency response 



can be determined from the above relationship. 


k. The frequency response was then used to filter the input (i.e., the 
normally distributed numbers) to the linear system. The filtering was done by 
utilizing the following equation: 


M 




V^+K 


K=-M 


where is the t n P u t to the system, a^ represents the filter factors 

obtained by transforming the frequency response into the time domain, and 

is the output which represents a random time history having the desired shaped 
power spectrum. 

5- Four time histories were generated in this manner. Power spectrums 
were calculated for each in order to insure that the proper filter factors had 
been obtained. 

6. Once it was determined that the calculated power spectrums were essen- 
tially the same as the desired shaped power spectrums, the , digital random time 
histories of step 5 "were analyzed to determine their instantaneous mean and 
amplitude distributions. 

The procedure described can be used equally well for digitally simulating other 
random time histories having arbitrarily shaped power spectra. 
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RANDOM NUMBER GENERATOR 


Random numbers were obtained with a fixed-point pseudo random number gen- 
erator developed by the National Bureau of Standards (ref. 4). Each generated 
random number Rjj was obtained from the previous random number R n _j_ by 

taking the last 11 digits of the product RqR^_i where Rq = 5^ and 

N = 1, 2, 3, . . . . Numbers were then selected at random from the generated 
R]y. Only the first 6 digits p of the randomly selected 11-digit number were 
used in this investigation. Approximately 160 000 random numbers were selected 
in this manner each having an equally likely chance of occurring (i.e., uniform 
probability distribution). The set of numbers obtained were all greater than 
or equal to zero but less than or equal to 999 999* 


TRANSFORMATION TO NORMAL DISTRIBUTION 


The random numbers were transformed into a normal distribution with mean 
equal to zero and unit variance by an approximate equation developed by Tukey. 
(See ref. 5*) This transformation was made in order to simulate a stationary, 
Gaussian random process. The transformation requires that the random numbers 
be between zero and one. Therefore, all numbers p were divided by 10^ and 
designated q. Tukey’ s transformation is 


X' 



( 1 ) 


where q is the modified random number and X' is the normal deviate. It was 
found in this investigation that when X' became greater than 2.4, there were 
significant departures from the normal distribution. Hence, it was necessary 
to use a corrected normal deviate X, as follows: 


X = X’ 


x = X’ + -|1-(0.13)(X’ - 2.4) 2 


(|X'| ^ 2.4) 

(I X’| > 2.4) 


( 2 ) 


It should be noted that these equations restrict the normal deviates to the 
range -5«73 ^ X ^ 5*73 which is no great handicap. 

A power spectrum was calculated by using equations (Bl) to (B3 ) and the 
first 40 000 normally distributed random numbers X. Due to storage limita- 
tions in the computer used, only 5000 numbers could be handled at one time. 
Therefore, 8 power spectra were calculated for each of the first 8 groups of 
5000 numbers generated. The 8 power spectra were found to be essentially flat 
and varied only slightly from each other, indicating that the sample size of 
5000 numbers was sufficiently large. An average power spectrum was obtained 
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from the 8 groups of numbers (white noise) and the remaining properties were 
calculated based on this average spectrum. 


FILTERING OF RANDOM NUMBERS 


The generated random numbers X, which when taken at discrete uniform 
intervals of time define a time history having a flat power spectrum, can be 
modified by numerical filtering techniques in order to change their amplitude 
response characteristics and thus change the power spectrum of the time history. 
The amplitude response characteristics can be changed by utilizing the input- 
output relation of power- spectral analysis, which states that the product of 
the input power spectrum 0i n (f ) and the square of the amplitude response 
. . ,2 



spectrum 


(sometimes called a transfer function) is equal to the output power 
0out ( ^ 





(3) 


The amplitude response 



can be determined from this equation since 


0i n (f) is the flat power spectrum obtained by calculation from the generated 
random numbers and 0 O ut( f ) is the specified or desired power spectrum. The 


amplitude response defines the changes that have to be made in the frequency 
domain in order to obtain the desired shaped power spectrum. These changes can 
be reflected in the time domain by taking the Fourier transform of the ampli- 
tude response. A time history comprised of discrete values Y± and having the 
desired shaped power spectrum can be calculated with the use of the following 
equation: 


*i 



K=-M 


a K^i+K 




where y i+K - 0 when i < M. The Fourier coefficients result from the 

Fourier transform of the amplitude response and the generated random numbers X 
are represented by y±±jf Details for determining the coefficients of the 

Fourier cosine series representation of the amplitude response are given in 
appendix A. The four amplitude response functions used in this investigation 
are shown in figure 1. The symbols show the shape of the response actually used 
to filter the random numbers whereas the solid curve shows the desired response. 
Twenty points were used to represent this response. These four amplitude 
response functions represent the concepts of bandwidth-limited white noise, 
atmospheric turbulence phenomena, single-degree-of -freedom system, and a modi- 
fied single-degree-of-freedom system (band pass), respectively. For brevity the 
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filtered time histories obtained by using these response functions are referred 
to as time histories A to D, respectively. Statistical samples showing the 
characteristically different features of the four time histories obtained in 
this manner are shown in figure 2. For clarity, the values of Y± have not 

been plotted but rather the curves faired through these values. The increment 
of time is At = l/2 sec between values of Yq. As a reference, 10 At is 
shown in figure 2. 

For the normally distributed numbers a power spectrum was calculated for 
each set of filtered random numbers using equations (Bl) to (B3) and the first 
k-0 000 numbers in each set. The power spectra were obtained by averaging the 
power spectra of 8 groups of 5000 numbers. Power- spectral properties were cal- 
culated based on the average power spectrum. This procedure was followed in 
order to determine whether the filtered time histories had power spectra 
equivalent to the specified or desired power spectra. The calculated power 
spectra were equivalent, within small tolerances, to the specified power 
spectra. 

In calculating these power spectra, the assumption was made that the fil- 
tered random numbers represented a sampling from a continuous time history 
y(t) at discrete uniform intervals of time At = l/2 sec which resulted in a 
discrete set of values Yi when t = i At. There is no loss of information 
from this sampling if the time history y(t) contains no frequencies greater 
than the Nyquist or folding frequency fp where 

" rk (5) 

The frequencies f, (2fp ± f), ( 4fp ± f), . . . cannot be distinguished in 
any frequency representation of y(t) which is determined from the values of 
yj _ . Thus, frequencies greater than fp will appear to be in the range 

0 < f ^ fp. It is said then that all the frequencies in y(t) have been 
folded into the range 0 ^ f ^ fp. This folding property follows from the 
relations 


sin 2jr(2mfp ± f )t = sin(&mti ± 2«ft) = ±sin 2rtft 


where 


i — 0, 1, 2, 3, . . . 

m = 1, 2, 3, 4, ... 

If frequencies higher than fp are present in the sampling, they will appear 

to contribute power or energy to the lower frequencies which will result in 
errors in the power spectrum at the lower frequencies. This situation was auto- 
matically eliminated by properly selecting the frequency range of the shaped 
output power spectrum. 
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POWER- SPECTRAL-DENSITY CHARACTERISTICS 



OF RANDOM TIME HISTORIES 


In using the power- spectral-density approach for analyzing the fluctua- 
tions of a random process it will be assumed that the process is stationary 
(i.e., statistical properties are invariant with time) and also Gaussian in 
nature. The power spectrum or power spectral density 0(f) is a frequency 
distribution function which describes the frequency content of the time varia- 
tion of a random disturbance y(t). For stationary processes, the power spec- 
trum 0(f) may be defined by the relationship which exists between 0(f) and 
the covariance or autocorrelation function R(t). This relationship is 
expressed as a Fourier cosine transform pair as follows (ref. 6): 


0(f) 


R(t) 



( 6 ) 


The covariance function, which is the mean value of the product y(t)y(t+T), 
gives a measure of the correlation between values of y(t) separated by a time 
interval t. Hence 


R(t) = y(t)y(t+T) = lim 

T-» OO 


1 _ 

2T 


f y(t)y(t+r) dt 


(7) 


For the special case when t = 0 

R(0) = y(t) 2 = f 0(f) df = ct 2 (8) 

u q y 

The function 0(f) may be regarded as the contribution of any frequency f to 
the mean square of y(t). The square root of the mean square value is known as 
the root mean square (rms) or standard deviation Gy of y(t). 

The derivatives of a Gaussian random process are required in order to 
determine* the statistics of such quantities as the number of times per unit 
time the disturbance crosses the axis y(t) =0, the number of maxima of y(t) 
per unit time, or the number of times per unit time that the disturbance 
exceeds a value of y(t) = y-^ where i = 1, 2, 3, . . . . The following rela- 
tionships are the ones developed by Rice (ref. 2) between the derivatives of 
y(t) and the power spectral density 0(f): 
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( 9 ) 


and 



( 10 ) 


The relationships involving these derivatives in obtaining zero crossings, 
peaks, and level crossings are as follows. 


The number of times per second that the zero axis is crossed with a posi- 
tive slope is 


r 

® 2jt tfy 


(H) 


The number of positive peaks per second is 


f _ i a y 


(12) 


The number of times per second that a value of y(t) = is exceeded is 


fy = f o ex P 


'll 

2 4 


( 13 ) 


(See ref. 6 for limits on y^.) 

Another relationship involving the derivatives of y(t) can be used to 
obtain the probability Pp(y ) that a peak will exceed a given value of 
y(t) = y^. The probability is expressed in terms of a standard variable z 
where z = y±/Oy 


The expression for the probability of obtaining a peak greater than y^ 
is (ref. 7) 


V y) - p i(i|) 



w 


" here P n(^) “ a P H(^) 

be exceeded; that is 


are the normal probabilities that 


— and — will 
% K 2 
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exp 


dz 


(15) 


where 


or 



jl _ r°° 

^ 2 rt J 7,/YL 


-Z 


2 


2 



(16) 


(17) 


The previous relations are valid only Tor a stationary random process 
which is Gaussian in nature. A more detailed discussion of this subject may 
be found in references 6 and 8. For data-processing purposes, the operations 
representing these expressions are more conveniently expressed in other forms. 
Appendix B gives the expressions which are in a form amenable to digital com- 
puting (eqs. (Bl) to (b6)). 


ANALYSIS OF RANDOM TIME HISTORIES 


An attempt is made in the present investigation to describe analytically 
the distributions of instantaneous means and amplitudes of several random time 
histories having different power- spectral properties and to develop relation- 
ships between these analytical expressions and the power- spectral properties 
of the random time histories used. This was done by analyzing four random time 
histories having different statistical and power- spectral properties which were 
generated with the aid of a digital computer. Some of the power- spectral prop- 
erties of these time histories have been calculated and are given in table I. 

The distributions of instantaneous means and amplitudes were obtained by actu- 
ally counting each mean and amplitude in the time histories. The number of 
occurrences of each of these values is listed in tables II to V. A discussion 
is presented of the distributions obtained by counting, the manner in which 
these distributions were described analytically, and the relationship between 
these distributions and the power spectra of the various time histories. 

The frequency distributions of the means for the four time histories inves- 
tigated are plotted in figure J. All four distributions appear to be normally 
distributed. This normality was checked by plotting the probability of 
exceeding a given mean value on normal probability paper for each of the time 
histories (fig. 4). As a first approximation the means can be considered to be 
normally distributed. 

The frequency distributions of the amplitudes f a about specified means 
are tabulated in tables VI to IX. Only those distributions which were 
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considered to have a sufficient sample size to be representative have been 
plotted in figures 5 to 8. The positive and negative amplitudes are approxi- 
mately symmetrical. Amplitudes were considered to be either positive or nega- 
tive depending on whether the slope of the line between successive peaks was 
positive or negative. The amplitude distributions are approximately symmetri- 
cal about the zero mean. 


An equation, similar to the one developed by Rice (ref. 2) for determining 
the peak distribution, which gave the best fit to the data was found to repre- 
sent the distributions of the amplitudes about any specified mean. The equa- 
tion represents the sum of a normal distribution and a Rayleigh distribution 
and is expressed mathematically as follows: 


fr 


“ % P N 


N r P r 


(18) 


where 

f c computed number of occurrences of amplitudes in range ^ y < 

N]\j number of amplitudes normally distributed 

Np number of amplitudes distributed according to Rayleigh distribution 

pjj no rmal probability which may be expressed by the following expression: 


% ~ 


^i+1 

I 

cr N \/at -y ± 


a- r 

J v . 


-a^ A 
exp da 


( 19 ) 


fr 


Rayleigh probability which may be expressed by the following 
expression: 


Fr - 


i+i 

4 


-p 2 


dp 


2d 


( 20 ) 


R 


The general form for equation (l8) is a modification of the peak prob- 
ability distribution equation which is the stun of a normal and a modified 
Rayleigh distribution. The coefficient op in the Rayleigh portion of the 
equation is the slope of the straight-line portion of the curve obtained from 
a plot of log of the cumulative frequency of amplitudes against the square of 
the amplitude for any specified mean. It was found that the coefficient ct r 

remained approximately constant, regardless of the mean value. The following 
relationship was developed, by using a trial and error procedure, in order to 
relate the coefficient dp to some of the power-spectral characteristics of a 

random time history: 


12 


( 21 ) 


a R = E 


+ (jjjj ~ 


K J 2 


+ f T 


where f Q , fp, Kq, and K 2 are derivable from the power spectrum of a random 
time history. 

The coefficient cr^ in the normal portion of the equation which gave the 
best fit to the data was found to be related to the coefficient or as follows: 

°N = 1 “ °r (22) 

The remaining coefficients in equation (l8), namely Nr- and Nr, were 

adjusted by a least-squares technique to give the best fit to the data 
(eqs. (C3) and (C4)). These coefficients might have some physical significance 
but for the purpose of the present paper they will be treated as being inde- 
pendent. The technique used is summarized in appendix C. Since the distribu- 
tions of amplitudes are approximately symmetrical, only the positive amplitudes 
were used to determine the coefficients. The number of positive amplitudes at 
a specified mean Nr was found to be related to the coefficients % and Nr 

by the following expression: 


% = n r + \ % 

which may be rewritten as follows: 

1 = % + 1 % 

Niji 2 Niji 

The ratio of Nr/^ was found to be a nonlinear function of the mean, being 

fairly constant for means close to zero and becoming progressively smaller for 
larger means. An attempt was made to predict the quantity % by determining 
the probability of occurrence of the means and multiplying it by the total num- 
ber of occurrences in the time history. Since it has already been established 
that the probability distribution of the means is approximately normal, all 
that is required is the standard deviation of the means. A relationship between 
the standard deviation of the means and the power- spectral characteristics of a 
random time history was found but the relationship was not sufficiently accurate 
to predict small standard deviations - that is, time histories C and D - and 
therefore not accurate enough to predict Np. No apparent relationship was 

found between the power- spectral characteristics of a random time history and 
the coefficients %j- and Nr. 

The coefficients derived to give the best fit to the observed frequencies 
for each of the time histories investigated are presented in table X. The 


(23a) 


(23b) 
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computed frequencies based on these coefficients are given in tables VI to IX, 
and the observed frequencies are also given for comparison* In addition, the 
observed frequencies (open symbols) and the computed frequencies (solid symbols) 
are plotted in figures 5 to 8. 

In the present paper, a strictly empirical approach was taken. An equa- 
tion was fitted to the data using a least-squares technique and two variables, 
namely, % and Nr# Possibly a better fit could be achieved by adjusting 

the four coefficients eft, cr R , N^, and Nr simultaneously with a least- 

squares procedure. However, a strictly analytical approach would be more 
desirable. It is most probable that an expression could be derived analyti- 
cal! y since the expression for amplitudes developed in this paper is quite 
similar to the expression for peaks developed by Rice (ref. 2). In addition, 
a definite relationship exists between the peaks and the means and amplitudes. 


CONCLUSIONS 


Four random time histories with significantly different statistical and 
power- spectral properties have been generated with the aid of a digital com- 
puter. The statistics of the means and amplitudes as well as the power- spectral 
characteristics have been obtained for each time history. The following con- 
clusions have been drawn from an analysis of the data obtained: 

1. The frequency distributions of the means are, in first approximation, 
normally distributed and symmetrical about a mean of zero. 

2. The frequency distributions of the positive (or negative) amplitudes 
for a specified mean can be described by the sum of a Rayleigh and a normal 
distribution. The positive and negative distributions are approximately sym- 
metrical. These distributions are also approximately symmetrical about a mean 
of zero. 

3- The standard deviations of both the normal and Rayleigh distributions 
representing the frequency distributions of the amplitudes are essentially 
constant over the entire range of mean values and can be approximated from the 
power- spectral characteristics of the time histories. 

4. The coefficients % and Nr in the general equation defining the 

distribution of amplitudes have been obtained empirically but no apparent rela- 
tionship between these coefficients and the power- spectral properties of the 
time histories has been found. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 3 , 1964. 
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APPENDIX A 


DETERMINATION OF COEFFICIENTS OF A FOURIER COSINE SERIES 
REPRESENTATION OF FREQUENCY-RESPONSE FUNCTION 


Consider the continuous periodic function of time, amplitude A, and fre- 
quency f, such that 


y(t) = A cos 2nft 

If this function of time is sampled at discrete time intervals At, the con- 
tinuous time history is replaced by a discrete set of values y-j_ equal to 
y(t) when t = i At and undefined in between. Thus, 

y(t) = y^ = A cos 2jtfi At (t = i At) 


When At = — i— . 

2f F J 


y ± = A cos ijt — 
f F 


(Al) 


The above time history can be modified to change its frequency character- 
istics (i.e., numerically filtered) as follows: 


*1 



K=-M 


a Ky i+ K 


(A2) 


vhere 


Y-^ filtered time history 

y. original time history 

i+K 


ag filter factors or coefficients 

M number of points used to approximate the amplitude response 


Equation (A2) represents, in numerical form, the passage of an input signal 
y(t) through some linear system which results in an output signal Y(t). Upon 
substituting y^ from equation (Al) for y 1+ ^- in equation (A2), 


Y^ = A/aQ cos tc 


f_ 

f F 




- K) + a^- 


cos tc — (i 

f F 
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Where a^- = a_g- 


/ 


Yi = A cos at — i 

fp 


a Q + 2 


V 


M 

I 

K=1 


cos at — K 


Yi = y ± a Q 


a n + 2 


M 

I 

K=1 


air cos at — K 

K y f 


(A3) 


The term in parentheses is in the form of a Fourier cosine series. It is 

necessary therefore to represent the amplitude response function F( — 

\ f F 

In the form of a Fourier cosine series in order to filter the generated time 
history. Therefore let 

_f_ = h 
f F H 


where h = 0, 1, 2, . . . H, then 


Ff|"l 

f F, 


H 

■) - F (|)- F h- a 0 + 2 X 

h=l 


"K 


cos at 


Kh 

H 


where 


®K 


-5 L 


H 


_ Kh 

F h cos at — dh 


Using the trapezoidal rule of numerical integration 


where 



1 

H 


cos at 




0 0 • • • 

1 

1 

1 

2 


1 6 



APPENDIX B 


I 


EXPRESSIONS USED FOR DIGITAL COMPUTING 

Given a time history y(t), digitized at discrete uniform time inter- 
vals At, and assuming that the origin of time occurs at one of these inter- 
vals, then 

t = i At 

where i = 0, 1, 2, . . . N and 

y(t) -* y(i At) = y ± 

The covariance or autocorrelation function shall be defined as a quan- 
tity Rp where 

N-P 

H P - „ + i . X Vwp (p - 0, 1, 2, . . . M - 60) (Bl) 

1=0 

Equation (Bl) is the numerical integration of equation ( 7 ). 

The power spectral density is the Fourier cosine transform of the covar- 
iance function Rp. For convenience it is obtained in two steps. The prelim- 
inary step gives estimates of the raw spectral density Ifo, and the final step 
gives estimates of the smoothed spectral density Estimates of spectral 

density are termed raw when they are obtained from the covariance function Rp 

by Fourier cosine series transformation and smoothed when hanned (operation of 
smoothing with weights 1/4-, l/2, l/b) from the raw estimates (ref. 9)* The 
smoothing operation partially accounts for the fact that a finite sample rather 
than an infinite sample was used when taking the Fourier transform. There are 
as many estimates as there are terms in Rp, that is, M. The raw spec- 

tral density estimates are given by the matrix equation: 
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where h is used to represent a frequency 

•p _ h r _ h 
11 M F 2M At 


and 


p = h = 0, 1, 2, . . . M = 60 

The smoothed spectral density estimates are given by the matrix equation 

[2 2 0 .. 

112 1 




0 12 1 




(B3) 


12 1 
0 2 2j 

The mean square of y(t) and its derivatives are defined as follows: 
°y = f>0 + fil + ?>2 + ' * * 0M-1 + \ 0m) 


Cb4) 


where Af^ = 


h 


2M At 2M At 


(aof ty 

PM Af.^l Z_j 



(B5) 



(at ) 4 j 

(2M At) 5 



(B6) 
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APPENDIX C 


LEAST-SQUARES TECHNIQUE FOR DETERMINING 
THE COEFFICIENTS % AND Njj 


The equation chosen to represent the frequency distribution of amplitudes 
about a specified mean is 

f c = % P N + N R P R (Cl) 

where f c is the computed number of occurrences of amplitudes in the range 
yj_ < y % y^+q* Pn or< l e r to facilitate computation of the number of occurrences 

of amplitudes in the range y > y^, fy^ is computed first . The desired value 
f c can then be computed from the following relation: 

f c - f y ± “ f y i+1 


Thus, 


f yi = % P N( y i) + N R P R(yi) 


(C2) 


where 


p N(yi) = 


i: 


lOO O 

-a 

exp da = 1 - 

2 o^ 2 a N \l& x 


r 1 exp da = 1 11 - 0(x)] 

■\]2rt J 0 2 ax 


where 0(x) = Error function 


i r“ -p 2 _y i 

p R(yi) = — p / P ex P — 3 d P = ex P — p 

0-R 2 J y± 2 a R 2 2o R 


The following approximation, obtained from reference 10, was used to facilitate 
computation of the error function 0(x) in the computer: 

0(x) = 1 - 


^1 + a^x + agx 2 + a^x^ + a^x^ j 
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APPENDIX C 


where 


x.SL 

2(Tn 

a± = 0.278393 
a 2 = O.23O389 
= 0.000972 
a^ = O.O78IO8 


The least-squares technique involves minimizing with respect to each of the 
undetermined coefficients the sum of the differences squared between the 
actual and calculated number of times a value y = y^ is exceeded - that is. 


Z { fh ' fy i) 2 ' 


Since the coefficients on and Op have been previously- 


determined to be constant it is only necessary to minimize with respect to Np 
and %. Thus, 



8 

8N N 



%%(yi) 



= 0 


= 0 


Solving these two simultaneous linear equations for Np and Np gives 


Np 


00 


00 


00 


00 


lM y i )] 2 I f * P R( y i) - Z M y i) F N( y i) Z 



00 00 

IW 8 lH4 2 

1=0 i=0 


Z P F( y i) P H( y i) 

i=0 


(C3) 


00 00 00 00 


Z f aM y i) Z P R 

«r. i=0 1=0 

:( y i) 

%(yi) - 

■ Y* faP N( y i 

i=0 

■) Z[ P *( y i)] 

i=0 

N “ 

00 

Z pr ( : 

_i=0 

^i) p n( 

y i) 

2 

00 00 

- Z[ p "(n )] 2 I[ p «(n 

i=0 i=0 

r 


(c4) 


20 


REFERENCES 


T .r 




4 


1. Schijve, J.: The Analysis of Random Load-Time Histories With Relation to 

Fatigue Tests and Life Calculations. Fatigue of Aircraft Structures, 

W. Barrois and E. L. Ripley, eds.. The Macmillan Co., 1963 , pp. 115-149- 

2. Rice, S. 0.: Mathematical Analysis of Random Noise. Bell System Tech. J. 

Pts. I and II, vol. XXIII, no. 3, July 1944, pp. 282-332. 

Pts. Ill and IV, vol. XXIV, no. 1, Jan. 1945, PP- 46-156. 

3- Naumann, Eugene C.; Hardrath, Herbert F. ; and Guthrie, David E.: Axial- 

Load Fatigue Tests of 2024-T3 and 7075-T6 Aluminum-Alloy Sheet Specimens 
Under Constant- and Variable-Amplitude Loads. NASA TN D-212, 1959- 

4. Taussky, Olga; and Todd, John: Generation and Testing of Pseudo-Random 

Numbers. Symposium on Monte Carlo Methods, John Wiley &. Sons, Inc., 
c.1956, PP* 15-28. 

5 . Tukey, John W. : The Practical Relationship Between the Common Transforma- 

tions of Percentages or Fractions and of Amounts. STRG Tech. Rept. 

No. 36 (Contract No. DA 36-034-0RD-2297), Dept. Math., Princeton Univ., 
June i 960 . 

6 . Bendat, Julius S.: Principles and Applications of Random Noise Theory. 

John Wiley & Sons, Inc., c . 1958 • 

7 . Huston, Wilber B.; and Skopinqki, T. H.: Probability and Frequency Char- 

acteristics of Some Flight Buffet Loads. NACA TN 3733, 1956. 

8 . Davenport, Wilbur B., Jr.; and Root, William L.: Random Signals and Noise. 

McGraw-Hill Book Co., Inc., 1958 . 

9 . Blackman, R. B.; and Tukey, J. W. : The Measurement of Power Spectra. 

Dover Publ., Inc., 1959- 

10. Hastings, Cecil, Jr.: Approximations for Digital Computers. Princeton 

Univ. Press, 1955- 


J 




21 


L 



TABLE I 


POWER-SPECTRAL- DENSITY CHARACTERISTICS OF THE 
FOUR RANDOM TIME HISTORIES 



TABLE II 


FREQUENCY OF OCCURRENCE OF IN STAM 1 AJTEOUS MEANS AND INSTANTANEOUS AMPLITUDES FOR TIME HISTORY A 


0.0 - 0.2 + 0.2 - 0.4 + 0.4 - 0.6 + 0.6 


-1.0 +1.0 -1.2 +1.2 -1.4 +1.4 -1.6 +1.6 -1.8 +1.8 -2.0 +2.0 -2.2 +2.2 -2.4 +2.4 -2.6 +2.6 Total 


53 42 28 27 15 20 4 

43 44 28 29 16 15 3 


44 56 22 25 22 15 3 

50 56 35 25 15 9 5 


60 69 29 16 13 13 7 

56 62 32 26 19 18 4 


69 57 26 26 i4 15 7 

60 78 24 42 l4 16 7 


61 76 36 32 12 17 8 

59 59 26 35 22 12 5 


65 70 26 24 8 20 4 

68 66 34 26 11 7 3 


47 53 23 22 9 7 4 

59 43 22 23 13 13 6 


4l 4o 16 15 

4l 42 15 15 


18 28 13 14 

24 30 8 11 


. 6,865 6,663 6,698 5,662 5,634 4,i46 4,291 2,810 2,931 1,726 1,767 980 1,031 465 458 227 234 85 80 24 20 8 10 2 1 
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FREQUENCY OF OCCURRENCE OF INSTANTANEOUS MEANS AND INSTANTANEOUS AMPLITUDES FOR TIME HISTORY B 



+ 1.1 

-1.5 

+1.3 

-1.5 

+1.5 

-1.7 

+1.7 

-1.9 

+1.9 

- 2.1 

+ 2.1 

-2.3 

+2.3 

-2.5 

+2.5 

-2.7 

+2.7 

-2.9 

+2.9 

-5.1 

+3.1 

-5.5 

+5.5 


187 15 k 159 122 172 116 125 


154 123 124 107 117 

130 11& 146 109 


104 

84 

60 

69 

50 

59 

23 

15 

6 

8 

11 

4 

4 

2 

1 

2 


82 

71 

64 

60 

56 

56 

21 

25 

13 

12 

12 

8 

3 

3 


89 

100 

48 

62 

39 

42 

29 

24 

9 

n 

6 

6 

3 


49 

88 

58 

59 
28 
37 
21 
13 

5 

10 

4 

4 

3 

7 


90 

74 

84 

44 

48 

29 

35 

19 

21 

11 

8 

2 

5 

l 

5 


103 

84 

72 

60 

48 

50 

28 

36 

21 

8 

6 

8 

1 

7 

2 

1 

1 


94 

89 

57 

70 

34 

39 

24 

43 

12 

24 

9 

7 

6 

6 

3 

1 

1 


76 

47 

73 

36 

47 

36 

28 

20 

13 

11 

6 


73 

72 

61 

45 

44 

27 

36 

22 

16 

5 

8 

5 

6 
5 
2 
1 
1 


52 

4o 

45 

22 

37 

16 
20 

9 

17 
6 
2 
2 
3 
3 

l 


57 

42 

57 

54 

29 

22 

27 

11 

9 

4 

10 

4 

2 

1 

l 

l 


39 

29 

27 

21 

17 

14 

5 

7 

10 

2 


33 24 


30 

4l 

22 

17 

16 

4 
7 

5 
4 
1 
2 
2 
1 
1 


16 18 


4 
10 

6 

5 

4 
1 
1 

5 


11 

9 

8 

7 

2 

5 

1 

3 


1,458 

1,125 

1,134 

765 

782 

513 

530 

326 

356 

186 

170 

78 

84 

52 

46 

18 

23 

8 

4 

3 

3 


Total 3,881 3,811 3,890 3,362 3,51a 2,874 2,860 2,133 2,138 1,497 1,599 1,074 1,086 607 637 362 389 202 242 107 122 57 57 20 17 


2 '36,561 



TABLE IV 



- 0.1 


+.3 

-.5 

+.5 

-.7 

+.7 

-9 

+.9 

- 1.1 

+ 1.1 

-1.3 

+1.3 

-1-5 

+1.5 

-1.7 

+1.7 

-1.9 

+1-9 

- 2.1 

+ 2.1 

-2.3 

+2.3 

-2.5 

+2.5 

-2.7 

+2.7 

-2.9 

+2.9 

-3.1 

+3-1 

-3-3 

+3-3 

-3*5 

+3-5 

-3-7 

+3-7 

-3-9 

+3-9 

-4.1 

+4.1 

-4.3 

+ 4.3 

-4.5 

+4.5 

I Total 


FRECSJESCI OF OCCURRENCE OF INSTANTANEOUS WEARS AND INSTANTANEOUS AMPUTUBES FOR TIME HISTOHI C 


0.0 

-0.2 

+0.2 

-0.4 

+0.4 

-0.6 

40.6 

-0.8 

+0.8 

-1.0 

+1.0 

-1.2 

+1.2 - 

Hoi 

324 

364 

232 

224 

144 

132 

49 

46 

17 

22 

6 

2 

375 

344 

355 

249 

219 

133 

l4l 

56 

55 

18 

13 

6 

5 

356 

332 

308 

192 

195 

89 

91 

51 

28 

6 

9 

2 


334 

321 

310 

173 

199 

73 

89 

33 

25 

8 

6 

2 

1 

467 

374 

425 

218 

229 

80 

90 

30 

28 

3 

6 

3 

1 

443 

426 

437 

240 

221 

86 

117 

35 

36 

5 

7 

1 

1 

608 

462 

492 

286 

291 

315 

129 

26 

30 

n 

10 


1 

585 

453 

501 

264 

294 

124 

124 

31 

32 

3 

2 

1 


690 

509 

403 

278 

290 

142 

116 

38 

34 

11 

9 

1 

1 

642 

526 

496 

•318 

285 

124 

122 

40 

31 

15 

12 



598 

516 

522 

298 

322 

1 16 

138 

37 

37 

5 

6 



599 

527 

540 

291 

317 

121 

121 

33 

43 

10 

8 


1 

576 

444 

468 

262 

278 

H5 

126 

37 

38 

5 

8 

1 


541 

429 

448 

270 

290 

129 

130 

4o 

4l 

6 

7 

2 

2 

438 

340 

380 

244 

226 

108 

115 

55 

38 

7 

4 

2 

1 

445 

365 

391 

214 

232 

102 

112 

28 

39 

7 

8 

2 

3 

312 

254 

295 

165 

106 

78 

77 

18 

22 

3 

6 

1 


325 

261 

295 

163 

172 

71 

85 

34 

18 

5 

3 



224 

227 

176 

118 

142 

53 

60 

18 

8 

3 

6 

1 


245 

213 

185 

121 

126 

62 

65 

l4 

27 

3 

5 

1 


170 

162 

139 

97 

87 

39 

36 

10 

10 

4 

3 


1 

175 

179 

127 

94 

91 

35 

26 

16 

9 

2 

4 



104 

98 

94 

57 

60 

16 

19 

7 

9 

1 

2 



104 

85 

no 

36 

49 

22 

24 

12 

12 


1 


1 

73 

50 

69 

36 

39 

23 

21 

6 

7 


1 



64 

59 

57 

33 

32 

23 

18 

5 

9 


2 



42 

36 

40 

23 

22 

6 

13 

5 

5 

1 




46 

33 

35 

21 

25 

« 

7 

2 

3 

1 




20 

16 

17 

12 

8 

6 

9 

2 



1 



23 

25 

19 

14 

12 

11 

7 

2 

2 





13 

15 

18 

7 

3 

2 








10 

7 

11 

8 

4 

2 

5 




1 



4 

3 

10 

1 

3 









10 

4 

8 

1 

3 


2 







2 

1 

6 


1 


1 







3 

1 

3 

1 










1 

1 

1 

1 





1 







2 

1 










1 

1 


1 

1 


1 







10,069 

8,425 

8,657 

5,04o 

5,178 

* 

(VI 

2,369 

730 

723 

160 

172 

32 
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TABLE VI 


PREDICTED AMD OBSERVED FREQUENCIES OF OCCURRENCE OF INSTANTANEOUS AMPLITUDES 
ABOUT A SPECIFIED INSTANTANEOUS MEAN FOR TIME HISTORY A 


Mean -1.4 -1.2 -1.0 -0.8 -0.6 -0.4 


- 0.2 


0.0 


+0 • 2 +0 • 4 +0 • 6 +0 • 8 


Amp. _ 

class marina f c f a f c f a f c f a f c f a f c f a f c f a f c f a f c f a f c f a f c f a f c 


+ 1.0 
i 


+1.2 +1.4 


f f 
•“•a x c 


f f 
•■•a x c 


f f f f 
x a x c x a x c 


0.1 

• 3 

• 5 
.7 

• 9 

l.l 

1.3 

1.5 


28 28 45 46 74 87 156 153 192 185 263 243 301 325 328 321 313 325 
35 26 50 51 97 91 lW 152 216 213 291 280 322 342 325 348 379 342 
32 29 56 60 107 104 182 171 256 257 330 335 440 395 44l 407 391 395 
24 33 60 69 129 119 193 195 309 294 405 384 450 451 483 466 463 451 
26 33 59 69 116 120 195 197 296 297 j 371 388 443 456 492 470 444 456 
34 30 68 63 87 109 174 178 262 269 380 351 399 412 4ll i426 413 412 
22 25 59 52 71 90 132 147 188 222 260 289 312 340 337 351 318 340 


15 19 ^l 39 64 68 110 in 142 168 207 219 278,258 


1.7 

8 

13 24 '28 

^0 

48 

74 

78 

91 n7 125 15V 177 

181 

166 186 

164 

181 

135 

154 121 

n8 

74 

78 

49 

48 30 

28 n 

13 

1.9 

8 

8} 13 

18 

35 

31 

30 

51 

65, 77 

89 100 104 

n8 

109 

121 

99 

n8 

76 

100 

70 

77 

35 

51 

17 

31 

16 

18 

7 

8 

2.1 

2 

5 

11 

11 

19 

19 

22 

31 

34 

46 

46 

6l 

51 

71 

60 

74 

63 

71 

kl 

61 

31 

46 

25 

31 

15 

19 

7 

11 

4 

5 

2.3 

0 

3 

4 

6 

8 

11 

15 

17 

18 

26 

2 k 

3^ 

21 

ko 

36 

k 2 

24 

^0 

24 

34 

22 

26 

12 

17 

8 

11 

2 

6 

0 

3 

2.5 

0 

2 

1 

3 

5 

6 

5 

9 

8 

lk 

13 

18 

14 

21 

5 

22 

17 

21 

7 

18 

8 

14 

5 

9 

3 

6 

2 

3 

1 

2 

2.7 

1 

1 

2 

2 

1 

3 

1 

5 

5 

7 

2 

9 

9 

11 

8 

11 

8 

11 

5 

9 

2 

7 

4 

5 

2 

5 

1 

2 



2.9 







2 

2 



2 

4 

2 

5 

4 

2 

3 

5 

1 

4 

2 

3 



1 

l 





3.1 











1 

2 

1 

2 

l 

1 

l 

2 







1 

l 






235 266 251 258 


230 215 176 183 152 153 75 87 44 46 29 28 
312 293 235 213 155 152 123 91 56 51 25 26 
380 347 253 257 183 171 n7 104 62 60 26 29 
356 386 286 294 207 195 no n9 78 69 42 33 
387 388 326 297 197 197 113 120 59 69 35 33 
367 351 269 269 162 178 99 109 66 63 26 30 
275 289 212 221 145 147 104 90 43 52 23 25 
202 219 139 168 no in 53 68 


f a - observed frequency of occurrence (i.e . , number of amplitudes counted in an interval) 
f c - predicted frequency of occurrence (eq. (Cl)) 



TABLE VII 


PREDICTED AND OBSERVED FREQUENCIES OF OCCURRENCE OF INSTANTANEOUS AMPLITUDES 
ABOUT A SPECIFIED INSTANTANEOUS MEAN FOR TIME HISTORY B 


Mean 

-1.4 

-I 

.2 

-1.0 

-0.8 

-0.6 

-0.4 

-0.2 

0 

.0 

+0.2 

+0.4 

+0.6 

+0.8 

+1.0 

+1.2 

+1.4 

AmpT^\^ 






























Hi 

class mark 

fa 

f c 

fa 

f c 

•^a 

f c 

•^a 

f c 

f a 

f c 

f a 

fc 

f a 

±c 


fc 

fa 

fc 


fc 

f a 

f c 

fa 

f c 

f a 

^C 

f a 

fc 


m 

0.1 

74 

89 

132 

136 

183 

184 

211 

236 

286 

282 

268 

344 

380 

363 

386 

382 

364 

363 

338 

344 

290 

282 

234 

236 

182 

184 

139 

136 

100 

89 

• 3 

63 

62 

86 

104 

128 

143 

166 

187 

212 

237 

251 

288 

330 

310 

313 

320 

295 

310 

294 

288 

213 

237 

176 

187 

l44 

143 

92 

104 

43 

62 

• 5 

45 

36 

77 

72 

107 

102 

163 

137 

211 

187 

223 

225 

264 

250 

249 

250 

291 

250 

244 

225 

219 

187 

189 

137 

115 

:102 

91 

72 

35 

36 

• 7 

43 

j 27 

71 

61 

96 

88 il47 

120 

186 

171 213 

205 

243 

231 

256 

228 >245 

231 :2o4 

1 

203 

182 

171 120 

120 

75 

88 

73 

61 

37 

27 

• 9 

32 24 

57 

56 

77 

81 

94 

111 

165 

159 |178 191 

223 

215 

192 

212 I89 

215 !l91 191 152 ;159 100 

in 

92 

81 

52 

56 

25 24 

1.1 

24 20 

39 

48 

52 

69 

76 

94 

116 136 122 162 1154 183 187 

180 1159 

: i 

183 172 |i 62 125 ;i36 


94 

57 

69 

33 

48 

19 20 

1.3 

.13 : 16 

27 

37 

45 

53 

73 

73 

84 104 1109 125 

1 

118 

i4i 

130 

139 

146 

l4l 

90 .125 

89 'io4 

61 

73 

57 

53 

4i 

37 

20 16 

1.5 

'10111 

17 

! 26 

37 

37 

47 

51 

60 

73 

88 

87 

71 

99 

84 

97 

100 

99 

84 

87 

70 

73 

44 1 51 

29 

37 

17 

26 

5 n 

1.7 

5 

7 

: 5 

; 16 

20 

24 

28 

32 

50 

47 

59 

56 

60 

63 

69 

62 

62 

63 ' 48 56 

39 1 47 

36 

32 

27 

24 

4 

16 1 9 

7 

1.9 

3 

4 

10 

10, 17 

l4 

13 

19 

36 

27 

37 

33 

36 

37 

39 

37 

42 

37 

35 

33 

43 j 27 

16 

19 

9 

14 

5 

10 

■ 5 

4 

2.1 

2 

| 

2 

1 6 

! 5 

! 2 

8 

6 

! 10 

8 

15 

13 

18 

25 

t 

| 20 

15 

20 

24 

j 20 

| 21 

18 

. 24 

1 l 

1 15 ! 8 

10 

10 

8 

1 

5 

■ 3 

2 

2.3 

1 0 

1 

l' 3 

3 

1 

4 

5 

5 

8 

7 

10 

9 

! 12 

: 10 

| 

8 

10 

11 

1 

10 

8 

9 

7 

7 

6 

5 

2 

4 

2 

3 

1 

1 

2.5 

0 

1 

1 

1 

! 1 

! 1 

2 

1 

2 

7 

3 

4 

* 

8 

5 

4 

5 

6 

5 

1 

: 5 

4 

1 

6 

3 

2 

2 

' 0 

! 2 

1 

1 



2.7 

; 1 

t 

0 


! 

1 

1 

0 

1 

1 

1 

7 

2 

3 

2 

2 

2 

0 

2 

; 5 

1 

l 2 

1 

ll l 

1 

! 1 

1 





2.9 

• 

1 

j 





1 

0 


| 

i 0 

1 



2 

1 

0 

1 





! 1 

0 







3.1 

| 

; 


_[ 



1 




1 

-l - . 

0 



1 

0 

1 

0 

! 



! 

; 





>1 




f a - observed frequency of occurrence (i.e., number of amplitudes counted in an interval) 
f c - predicted frequency of occurrence (eq. (Cl)) 






TABLE IX 


PREDICTED AND OBSERVED FREQUENCIES OF OCCURRENCE OF 
INSTANTANEOUS AMPLITUDES ABOUT A SPECIFIED 
INSTANTANEOUS MEAN FOR TIME HISTORY D 


Mean 

-0.6 

- 

0.4 

- 

0.2 

0.0 

+0.2 

+0.4 

+0.6 

class marl 

•fa 

fc 

fa 

f c 

fa 


fa 



fc 

fa 

fc 

fa 

fc 

0.1 

73 

94 

231 

194 

395 

355 

486 

485 

355 

355 

194 

194 

94 

94 

•3 

48 

25 

149 

121 

340 

309 

491 

427 

341 

309 

154 

121 

46 

25 

• 5 

29 

38 

176 

184 

406 

470 

650 

650 

444 

470 

163 

184 

4o 

38 

• 7 

44 

k 6 

179 

225 

573 

576 

768 

795 

556 

576 

207 

225 

37 

46 

• 9 

43 

50 

227 

242 

647 

618 

903 

854 

649 

618 

236 
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f a - observed frequency of occurrence (i.e., number of 
amplitudes counted in an interval) 

f c - predicted frequency of occurrence (eq. (Cl) ) 
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Figure 3*- Statistical distributions of instantaneous mean values for the four time 

histories. 
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Figure 5.- Statistical distributions of instantaneous amplitudes with respect to a spec- 
ified instantaneous mean value for time history A. Open symbols represent actual 
values; solid symbols represent computed values. 
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Figure 5.- Concluded. 
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Figure 6.- Statistical distributions of instantaneous amplitudes with respect to a spec- 
ified instantaneous mean value for time history B. Open symbols represent actual 
values; solid symbols represent computed values. 
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Figure 6.- Concluded. 
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Figure 7.- Statistical distributions of instantaneous amplitudes -with respect to a spec- 
ified instantaneous mean value for time history C. Open symbols represent actual 
values; solid symbols represent computed values. 
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Figure 8.- Statistical distributions of instantaneous amplitudes with respect to a spec- 
ified instantaneous mean value for time history D. Open symbols represent actual 
values; solid symbols represent computed values. 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 
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